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Quantum cosmology of the multi-field scalar 
matter: some exact solutions 
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Abstract 

We consider the gravity interacting with matter scalar fields and 
quantized in the minisuperspace approach in which the wave func¬ 
tional is described by the Wheeler-DeWitt equations (WdW). Assum¬ 
ing the domination of the homogeneous and isotropic geometry the 
leading contributions to the wave functional in the approximation of 
the minisuperspace with Priedmann-Robertson-Walker metric (FRW) 
and spatially uniform scalar fields are considered. The model of several 
scalar fields with exponential potentials and kinetic terms admitting 
such a special mixing that ultimately it is possible to separate the 
variables in the WdW equation and to find its exact solution in terms 
of the special functions is proposed. The semiclassical approximation 
is thoroughly investigated and the boundary conditions permitting the 
physical solution selection for classical cosmologies are chosen. 


1 Introduction 

The exact Einstein eqnation solntions in the presence of the matter play 
important role in the cosmology. Last decades the relativistic models with 
the scalar matter that are applicable to describe the inflationary expansion in 
the early stages of the evolntion n.H and to describe the cosmic acceleration 
[3],[1] at later stages (the quintessence and the dark matter models [5]-[9]) 
have been quite popular. At the same time there are only few known exactly 
solvable cosmological models with scalar matter. One such a model is of a flat 

*V.A. Fock Department of Theoretical Physics, Saint-Petersburg State University, ul. 
Ulianovskaya, 198504 St. Petersburg, Russia 


1 



Friedmann universe filled with a scalar field with the exponential potential. 
Various solutions for one scalar held were studied in Hi-lii. Recently the 
general Einstein equation solution for this model were found |19] |33]. 

We mentioned classical solutions for the universe geometry evolution. 
However to understand in the details the metric singularity in the moment 
of the Big Band one has to know the quantum properties of the gravity 
interacting with the quantum matter helds. Without a doubt any exact so¬ 
lution for the quantum univers is invaluable even in the restricted sense - 
with few degrees of freedom quantized. To ensure unambigous classical limit 
the ADM quantization in the minisuperspace is commonly used practical 
essence of which is given by the Wheeler-DeWitt equations (see [31] and 
references therein) rehecting the constraints due to the symmetry under dif- 
feomorphisms. Generally speaking these functional equations are very hard 
to solve. However to some extent one can assume that for approximately 
flat universes the leading quantum effects are determined by the spacially 
uniform and isotropic fluctuations of the gravity and matter whereas the 
inhomogeneities and anisotropies can be succesfully analyzed in the approx¬ 
imation of the weak quantum fluctuations. In this case we come to the 
minisuperspace reduction with gravity and matter helds depending only on 
time. As a matter of fact such a model is a quantum mechanics with hnite 
number of the degrees of freedom. 

In this work we assume such a domination of the spatially uniform isotropic 
geometry and study the leading contribution to the wave functional in the 
FRW metric minisuperspace approximation with auxiliary lapse held that 
provides the time reparametrization invariance and scalar helds depending 
only on time. The model of several scalar helds with kinetic terms admitting 
mixing is proposed. This mixing is chosen in such a way that ultimately it is 
possible to separate the variables in the WdW equation and to hnd the exact 
solutions in the terms of special functions. In the section 2 we study the clas¬ 
sical Einstein-Friedmann euqations, solve them and thus provide the basis for 
consistent quantization of the WdW equation in the minisuperspace approx¬ 
imation (section 3). Exact quantum solution of this equation are obtained 
in the section 4. In the same section the boundary conditions for selection of 
the physical quantum solutions that admit reasonable classical cosmologies 
are introduced. Semiclassical (WKB) approximation is thoroughly investi¬ 
gated in the section 5. Future plans and prospects of this class of models are 
briefly discussed in the conclusion. 
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2 Model of several scalar fields with expo¬ 
nential potentials 


Let us consider the following model with several scalar helds minimally cou¬ 
pled to the gravity, 

^ ( 1 ) 

ab a 

where index a goes from 1 to n. 

Let us restrict our consideration to the minisuperspace corresponding to 
Friedmann-Robertson-Walker metrics and spatially uniform helds, 

ds'^ = (pa = ( 2 ) 


The corresponding minisuperspace action takes the form. 



It is convenient to search for the classical solutions in the Hamilton for¬ 
mulation. The canonical momenta are. 


= = Pn = 0. (4) 

b 

As always in the theory of gravity, the Hamiltonian becomes purely a con¬ 
straint with a Lagrangian multiplier N, 


H = Ne-^P 




abPaPb 




^Qp-\-Xa(pa 


ab 


(5) 


To eliminate the dependence of the exponent term on the different degrees 
one should to perform the following canonical transformation. 


\(Pa,PaJ \Xa,T(aJ 


6 

^ =Pp- y^^Pa, Xa 


Aa4a + 6p, 


. Pai 
'^a 


( 6 ) 

(7) 
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then the Hamiltonian takes the form, 


H = Ne-^P 


^2 2 
- UJ — K U 

12 


ab 


A„A 


+ ^)a 6 - 3 K^ 




E''«' 


Xa 


( 8 ) 


Note that new canonical momentum u is conserved on the constraint 
surface, i.e. 

{oj,H} = -ZH ( 9 ) 

Generally speaking this model is not exactly solvable in the case of several 
fields. Let us assume the special form of the kinetic term matrix. 


{M-^)ab = 


Dg^ab + 

Aa Ab 


( 10 ) 


such that the kinetic energy of the Xa fields becomes diagonalized. 


H = Ne 


-3p 


^2 2 
- U — K U 

12 




( 11 ) 


As will be shown below this permits separation of variables in the classical 
as well as in the quantum case. 

The classical solutions can easily be found by solving Hamiltonian equa¬ 
tions of motion. Let’s choose the gauge. 


N = e^P. 


(12) 

The equations of motion on the constraint surface take the form. 


Xa = {Xa,H] ~ -K^U + DaTTa, 

TTa = {vTa, iL} ~ -146^“, 

(13) 

P = {P: H) ~ -yCU - 

a 

W = {Xa,H} ~ 0, 

(14) 


Taking the derivative with respect to time of the first equation we get. 


Xa 


-DaVaC^^ 


( 15 ) 
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from which the following solution is obtained, 


gXa 


•iPI 1 

DaVa COsh^(P„t - Qa) ’ 


( 16 ) 


where Pa,Qa - integration constants. Substituting this result into other equa¬ 
tions we obtain general solution, 

TTa = -^tanh{Pat - Qa) + P = Po + Pa, ( 17 ) 

a 

pa = {Pat- Q a), (18) 

Pl = ^{Qn^ + Dan-^ + CaDa), E^7a = 0. (19) 


where some integration constants happen to be hxed thanks to the constraint 
equation H = t). 

Note that for solutions to be real requires. 



DaVa 


> 0 , 


( 20 ) 


which constrains the classically allowed values of Ca- For DaVa < 0 this 

~ 2 

condition demands P^ = —Pa <0, 


2P2 1 

6^“ = _2_^ ('21) 

\DaVa \ {Pat - Qa) ’ 

4p2 ^2 _ _ 

Pa = +-f^^t+ —\nCOS^{Pat- Qa). (22) 

For Pat — Qa = ±f the cosine vanishes and both Xa and p tend to infinity. 
Hence in this case the solution is to be considered on the finite interval of t. 
Note that in the case of > 0,DaVa >0 the classically allowed region of Xa 
is bounded from above while for P^, Pa 14 < 0 it is bounded from below. 

The DaVa < 0 case may correspond to the phantom fields (quintoms [35] ) 
or negative potentials unbounded from below. The question of the stability 
for such cosmological systems lays beyond the scope of this paper and will be 
left for the further analysis. Nevertheless it should be noted that according to 
(ITOj) sufficiently small negative value Da may still correspond to the positive- 
dehnite kinetic term matrix Mab- 
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Let us consider the case of one field with the different time variable anal¬ 
ogous to the one considered in [2I],|32]- 


dr 


dt DV 


P X 
62, 


For > 0,Fy > 0 we get, 

T = Tq + arctan tanh 


Pt-Q 


gXa _ 


Pa = 


4^ 

Da 


uj{Q + 2 arctanhtan(r — tq)) -|- 


2P2 1 

a _ _ 

DaVa cos2(r ■ 
— In cos2(r 



(23) 

T-o)’ 

(24) 

ro), 

(25) 


Note that original time variable t corresponds only to the finite interval of 
r. However at the ends of these interval the metric becomes singular = 0 
and beyond it becomes complex. 

For < 0,FH < 0, 


T = Tq + arctanh tan 

4Fa . , 

pa = 


Pt-Q 


eXa — 


2P2 


Da 


\DaVa \ cosh2(r - To)’ 
^2 

-In cosh 2( 

Da ^ 


T-o)’ 

(26) 

■ T-o), 

(27) 


he. the hnite interval of t in this case corresponds to the whole real axis of 


T. 

Cosmic time for which the metrics takes the form, 

ds^ = dtl^^^-e^Pd^, tcosm= [ dtN. (28) 

Jto 

Let us consider the case with all Da > 0,14 > 0 assuming that all helds 
pass the turning point approximately at the same time i.e. Qa — 0. For large 
values of |f| the metric factor p behaves as a linear function, 


P ~ do- 




K 


2VQDa 


a/ 6^2 + DaU 1 CaDa u\t\ (29) 


Let a; < 0. Then the cosmic time between the Big Band at f ^ —oo and the 
turning point at t = 0 is hnite. On the other hand the cosmic time between 
the turning point and the Big Crunch at f ^ -|-cxo is hnite only if, 


< 


-|- DaU ^ -|- CaDa 


1 + E 


6 k 2 

a Da 


(30) 
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To consider the solution with a; > 0 it is sufficient to reflect the time t ^ —t 
in the solution with negative u. 

If for at least one field Da < 0, 14 > 0 then the finite interval of t for 
which the solution is defined corresponds to the infinite cosmic time between 
the Big Bang and the Big Crunch. 


3 Wheeler-DeWitt equation 


The Wheeler-DeWitt equation for the minisuperspace under consideration is 
obtained by canonical quantization for the Hamilton function l|5|). 


Ne 


-3p 


^)abdadb + 

ab a 


,6p+Aa0a 


= 0 . 


(31) 

Factor Ne will not play the role in finding the solutions and thus it will 
be ommited. 

Canonical transformation (I6]),(17|) in the quantum case becomes simply 
the coordinate transformation and the corresponding transformation of the 
partial derivatives after which the equation takes the form, 


+ h^K^dp^da - ^)ab - ?>K^')dadl 


^n.^b 


ab 




'I’(ft{Xa}) = 0. 


(32) 


Because now there is no p in the potential, it is convenient to perform 
the Fourier transform. 


*(p. {Xo}) = J 


(33) 


^n.^b 


+ihK‘^UJ^ da - - Su'^^dadi 


ab 




5 X 0 


= 0 


(34) 
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In the case of the special form of the kinetic term matrix 


{M-^)ab = 


DgSab + 


( 35 ) 


mixed terms are canceled and partial solutions can be obtained by separation 
of variables, 

{Xa}) = n Aa), (36) 

a 

with H'a depending only on one coordinate Xa- 
The corresponding equations take the form, 


^ 2/ -1 

--2^ <" 


+ C„) + ihK‘ud^ - tt‘^81 + Ue*" 


^a(c<^,Xa) = 0, (37) 


where Ca are arbitrary constants under condition '^Ca = 0. 

a 


4 Exact solutions of WdW equation 


Let us consider now the exact solutions of the equations ([37j). 
variable change, 

kP'u} 

we obtain the modihed Bessel equation. 



After the 
(38) 


+ zd^ + pI 


f{u,z) = 0, 


4P 

p = ° = 

“ HD, ’ “ 24 


2 fa 2 


{6k + DaU + CaDa)) 


with Pa expressed similarly to the classical solution ffT6|) . ffT^ . 
Therefore we obtain the exact solution, 

Xa) = A^+)Ti+)(a;, Xa) + A^“)TW(a;, Xa), 


/ \ i K u 

= Qn Da r- 


2 2Va 


aXa/S 


fiV -Da 


(39) 

(40) 

(41) 

(42) 

(43) 















where take the following form in terms of modihed Bessel functions, 


L.{z) = (44) 
/sm(^ 7 r) V / 

Let us consider the case > 0, DaVa > 0. Consider the limit Xa —t —cxo 
{z —)■ 0+) which also corresponds to p —)■ —cxo for fixed (pa- Using the 
asymptotics of the modified Bessel functions with purely imaginary order 
near zero [52] we get that both particular solutions behave like plane waves, 


v]>(+) ^ _ 
a 


TT 


Ua Sinh UaTT 


Qli Da 


TT 

i^a sinh i^„7r 


Da 


sm I - 5a ) + O ( 
TZaXa 


COS 


A 1 1 ^ 


Sa]+Oi^e^' 

r(i + iv) 


(45) 

(46) 

(47) 


Let us consider now the limit Xa ^ +oo (z —>■ +oo) which also corre¬ 
sponds to p —-|-oo for hxed (pa- Using the leading asymptotics of Bessel 
functions Inu{z) and Knu{z) at large values of argument | 2 ;| ^ oo (according 
to 8.451 in PH) we obtain. 


\[/(±) 

a 


2 



V4 

e 4 exp 


2 fWa ^ 
^ hy Da""^ 


e^' 


-Xa 


(48) 


The selection of the solution may be done based on the condition of square 
integrability by p from —oo to -|-oo to ensure the Tip operator hermiticity 
and the self-consistency of the Dira or BRST/BFV quantization scheme [38] . 
Because this condition has to be satisfied for any value of {(pa} it yields the 
decreasing of the solution at large values of any Xa- Thus for the case of 
> O^DaVa > 0 , 


^{pAXa}) 





(49) 


Let us consider the case of = —P^,VaDa < 0. Then it is convenient to 
write the solution in the form, 

= AW¥X(^,Xa) + (50) 
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/_i_\ i K (jj .. 

X]>(±) ^ gR Ua J. 


±i>a 


2K 




pXa/2 


where Va = |i^a|- Now the limit Xa —t +oo {ya = | 
corresponds to the oscillations, 


2Va 


Da 


(51) 

gXa/2 ^ +00) 


^(±) 


TT V2Zln 


\ V4 

/2 /■ 

2K 

) e 4 

cos — A/ 

/ 

VhV 

L>a 


3 X 11/2 


TT 




sin 



2K 

Uv 

Da 


3X11/2 


vr _ 

2 “'“ 


-0+O(e-i-) 


(52) 

(63) 


whereas in the limit Xa ~oo {y —)■ 0+) (using the Bessel function asymp¬ 
totics from 8.440 in [2Z]), 





± 


T^aXa \ 
2 ) 


(54) 


Square integrability by p demands the selection of the decreasing solution, 

i.e. 


^(P,{Xa}) 


= j do; y4(a;)e'i‘^^ e'^' 


~Xa 


Jr, 



(55) 


Note that in both cases the selection was done according to the behavious in 
the classically forbidden under-barrier regieon. 


5 WKB asymptotics 

Let us consider WKB asymptotics of the Wheeler-DeWitt equation ([32D, 

{Xa}) = ^/^(p, (56) 

Hamilton-Jacobi equation takes the form, 

- K\d,S) J2{9rS) 

a 

+ - ^'^")idaS){dbS) + Y,Vae^‘^ = 0. (57) 

ab a 
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Similarly to the case of the classical and quantum equations the cancela¬ 
tion of the mixed terms permits the separation of variables, 

'5(p,{Xa}) = a;p +(58) 

a 


that corresponds to the separation of variables in the quantum equation, 

^(P,{Xa}) = ~ (59) 

a a 


Because u is conserved in the classical case these quasiclassical solutions 
are the ones that correspond to the correct classical limit. 

The Hamilton-Jacobi equation on Sa, 


- ^ +Ca)- K^uidaSa) + ^{daSaf + ^6^“ = 0, (60) 

has exact solution, 

s .=- xi“’)+ ~ ’ 

where €„ = ±1 and -|- DaU ^ -|- CaDa) and the constant is 

( 0 ) ‘2P'^ 

chosen for solution to vanish at the turning point e^“ = ,. ° . 

Treating the action fl^TD as the classical Hamilton-Jacobi action one can 
obtain the classical trajectories by differentiating it with respect to parame¬ 
ters uj and Ca- Then demanding that. 


dS 

—— = const, 

OOJ 


dS 


= const, 

EC'a=0 


(62) 


should yield the equations on the trajectory. Substituting the solutions 
flTGll . ffT^ and differentiating with respect to time one obtains. 


- (^) (l - sgn(Pat - Qa)^ = 0, (63) 

-4ea sgn(Pat - Qa) + 4ei sgn(Pit - Qi) = 0. (64) 


Thus the sign of is chosen depending on the direction of motion towards 
the turning point. 
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(65) 


The leading order of the prefactor ?/)„ happens to be equal, 


= (l 


VaDa 

2P2 


-1/4 


0{h). 


In the turning point WKB approximation ceases to work however one can 
match the solutions in the classically forbidden and allowed regions approx¬ 
imating the potential with the linear function V ~ VaXa and the solu¬ 


tions with Airy function 
Xa +00 it yields, 


For 


p2 

'D^ 


> 0 and the solution decreasing at 




^ ~ Xa ) ^ Tx/f — arctauh \/P -|- — y(66) 

4 a ^ 4 / 

—F — arctan \/—Pj ^ ,(67) 


B ijP 

— e f^Da 




a 

(0)^ / 4P„ 

exp 


(-P)i ^ V 

where B is an arbitrary constant and, 

VaDa 


P = 1 - 


Xa 


2P2 


( 68 ) 


Comparing this result with the asymptotics for Bessel function at large values 
of order HD] we obtain that it corresponds to the correct asymptotics of 
in m if, 


1 hnDa f.K (0) 


27rPa 

hDa 


(69) 

p2 

In the case of < 0 we have to choose the solution decreasing at 

FJa V a 

Xa -OO, 


g>(Xa>xi°D 


2P —v*-°b • 

— QhDa '■X“ Xa ) 


[-FY 


^(Xa<xi ^ y gCoa^Aa xi^lg^p 
Pi 


(^V-F - arctan V-F + 0, (70) 
^\/P — arctanh \/pj | , (71) 


hDa 


The following value of B corresponds to the asymptotics of in fl5T|) 
(using 8.452 in [37]). 


B = 


hDa 


2 V 27rP 


exp I z 


K 


hD, 


-^Xa 


( 0 ) 


(72) 
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6 Conclusion: prospects 

We have proposed the model of several scalar helds interacting with gravity 
that admits exact classical solutions thanks to the separation of variables 
in the special gauge (for special time variable choice). As well analytical 
solutions can be found after its quantization in the ADM approach for WdW 
equation. 

The quasiclassical correspondence of the solutions is investigated for the 
further definition of the physical time. In this paper we considered the special 
class of scalar field interactions both with positive and negative kinetic terms 
and potentials but only admitting real classical trajectories. Though all the 
classical solutions are found the question of stability for partially phantom 
systems (quintoms [35]) and negative potentials unbounded from below is 
left for further analysis. 

In the future we hope to investigate the classical solutions for PT sym¬ 
metric complex potentials [H] and compare the cosmological trajectories in 
these potentials with the trajectories in the systems with phantom scalar 
fields (quintoms [35]b 
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